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We numerically study the momentum distribution of one-dimensional Bose and 
Fermi systems with long-range interaction g/r'^ for the "special" values g = —^,0,4, 
singled out by random matrix theory. The critical exponents are shown to be inde- 
pendent of density and in excellent agreement with estimates obtained from c = 1 
conformal finite-size scaling analysis. 
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I. INTRODUCTION 



In one-dimensional (ID) quantum systems it has been argued that the existence of a 
critical point at zero temperature implies asymptotic conformal invariance at large distances 
p!9|,|3|l . This invariance powerfully constrains the possible behavior of a given critical theory, 



parametrizing universality classes by the conformal anomaly c, the value of the central charge 
of the underlying Virasoro algebra 0. For c < 1, it has been shown by Friedan, Qui and 
Shenker that unitarity (reflection positivity) restricts those values of c to be quantized [Q. 
The critical exponents are given by the Kac-formula and the correlation functions can be 
determined. Realizations of this scheme are the Ising (c = |), tricritical Ising (c = 
three-state Potts (c = |), tricritical Potts (c = |), and other models PJ^l- addition, a 
complete classification can also be given for c > 1 if a continuous symmetry is assumed; 
e.g. for f/(l) we get the Gaussian model (c = 1), describing a variety of models such as the 
Coulomb gas, XY-model, s = | antiferromagnet and more 

The basic assumptions leading to conformal invariance of a model may be summarized 
as scale invariance, rotational and translational invariance and short-range interactions. If 
any of the above properties do not hold, the conformal symmetry will in general not be 
obeyed either. In addition, it should be noted, that there is no rigorous proof of conformal 
invariance at a critical point, just as for scale invariance. Instead conformal invariance 
should be regarded as a principle, whose consequences can then be compared with exact and 
numerical results on specific models P]. In the case of long-range interactions, an explicit 
counterexample may be given: Thouless showed that certain long-ranged ID classical 
systems may have a critical point. Now, since any general coordinate transformation in one 
dimension is conformal, this would imply that all correlation functions are constant, which 
is certainly not true. 

Recently, however, a number of groups have independently applied conformal arguments 
to long-range models such as the periodic g/r'^ models |T^-|15|,|18[ . They were able to compute 
the conformal anomaly and the conformal scaling dimensions, as well as various correlation 
exponents for Bose and Fermi statistics as functions of the interaction strength g. Their 
predictions agree with the exact result for the case of g = 0, that is, free fermions and hard 
core bosons. For less trivial interaction strength, no independent verification of these results 
exists so far. 

The periodic g/r"^ model has been solved by means of the Bethe Ansatz by one of 
us several years ago |Q. As it turns out, for the three "special" values g = — 1,0,4 of 



the interaction strength, the ground state wave function of this model can be viewed as a 
distribution function for the eigenvalues of a random matrix from an appropriate ensemble. 
Application of theorems from random matrix theory then allows the derivation of explicit 
formulae for the ground state correlation functions p3[ . 

In this paper, we will present a high-precision study of the critical exponents for the 
momentum distribution, which we will then compare to the results obtained by the conformal 
conjecture. In section |l[ we will outline the arguments of conformal invariance leading to 
the predictions for the critical exponents. Section |T| introduces the g/r^ model and its 
ground state correlations in more detail. We present the results of our numerical study in 
section where in addition we present an unusual scaling law. 
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II. DERIVATION OF CRITICAL EXPONENTS 



Finite-size scaling analysis of classical two-dimensional statistical systems at criticality 
may be used to obtain the value of the conformal anomaly c and the scaling dimensions 
Xn- For short-range interactions that give rise to a linear dispersion relation, i.e. gapless 
excitations, the following formulae hold [Q: for a strip of width L — oo and periodic 
boundary conditions, 

Eo^eoL-^^, (1) 
o L 

En-Eor^ 2ir^. (2) 

Here En is the eigenvalue of the Hamiltonian; Eq is the ground state energy; eo is the ground 
state energy density in the thermodynamic limit, and v is the velocity of the elementary 
excitations. Since the partition function of a classical system of finite width and periodic 
boundary conditions may be viewed as the Feynman path integral for an infinitely long 
quantum chain at finite temperature, these results apply to ID quantum systems as well 



We now wish to evaluate the scaling dimensions for the ID Bose and Fermi systems 
interacting via the long-range interaction of the periodic g/r"^ model. This model has been 
solved by the asymptotic Bethe Ansatz in p2[, and it has been shown that the energy 
spectrum may be written as 



M 



=1 L 



2n ttA , 

-/, + — (2J-M-1) 



(3) 



where {/,} denotes the set of ordered (free boson) quantum numbers of a particular state; 
M is the number of particles; L is the length of the ring and 2A(A — 1) = g. The ground state 
energy Eq corresponds to the state with quantum numbers Ij = for all j = 1, . . . , M. In 
addition, the model has a sound type excitation spectrum, and for A = |, 1, 2 all correlation 
functions may be calculated. We will elaborate on this last fact in more detail in the next 
section. 

Let us now calculate the finite-size corrections to the excited states. We start with the 
Bose system first. An excitation corresponding to a momentum transfer P = is given by 
changing the number of particles M to M + AM. By (§), the excitation energy is given as 

^2 \ 2 ( M+AM M 

E{M + AM) - E{M) = -— { Yl [2j-(M + AM)-l]^-^[2j-M-l]^ 

-TT^A^d^AM+VA^dAM^ 

= ^AM+^^AM' (4) 

Here we used the formulae for the chemical potential = ir'^X^if and the sound velocity 
V = 27rXd, as obtained in 1221 • 



Let us define the Fermi momentum as kp = nd. Then for Fermi and Bose systems 
and fixed particle number M, taking AD particles from the "left" Fermi momentum to the 
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"right" Fermi momentum — or equivalently, by translating the filled Fermi sea in momentum 
space — results in another type of low-lying excitation, i.e. 



M 



E{M; AD) -E{M;0) = ^ 



|aD + ^(2j-M-1) 



1 2 



ttA 



(2j - M - 1; 



2ttv 1 , 9 



(5) 



The momentum of this excitation is P = 2kFAD. 

For fermions, the same reasoning applies. However, the quantum numbers Ij are re- 
stricted to integers now and we have to follow the selection rule 



AD 



AM 



mod 1). 



(6) 



By combining (H) and (0), we obtain the scaling dimensions for the above excitations as 



xiAM; AD) = ^AM^ + j^D^- 



(7) 



This expression is valid for both Bose and Fermi statistics modulo the selection rule @. 

We have not calculated the effect of particle-hole excitations on x. Although we do 
not need the resulting additional terms for the calculation of the long-distance behavior of 
the correlations, we remark, however, that in order to correctly read off the right and left 
conformal weights that characterize the underlying Virasoro algebra, they cannot be ignored 

laniii- 

Before we continue with our derivation of the critical exponents of the momentum dis- 
tribution, let us pause and first calculate the critical anomaly c for the present model. From 
CH) it follows that 



E^ = E, 



{0} = (^oL - -Ti'^X^d/L. 



(8) 



Using the above expression for the sound velocity f , we may rewrite the L correction term 
as — A. Comparison with (Q) yields c = A. This is a peculiar result in as much as the 
critical anomaly varies continuously with the interaction strength A. We could instead use 



the low-temperature expansion of the free energy as given in and the formula for the 
critical anomaly as in giving 



F{T) ~ F{T = 0) - 



(9) 



and therefore c = 1. Thus the two expressions for c are in contradiction. We will discuss 
this point in more detail in section 0. 

Let us finally investigate the long-distance behavior of the correlation functions. Let 
the pair (AM; AD) describe an excitation where the particle number is changed by AM, 
and AD particles are transported round the ring. If we then denote by cf)^ (0) the creation 
(annihilation) operator of this excitation, the corresponding leading term of the equal-time 
correlator has the form 
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(0''^(r)0(O)) oc cos(Pr) r 



-2x{AM;AD) 



where P is the associated momentum carried by the excitation Therefore the correlation 
of the Bose field operator has the form 

(0UO0b(O))a cx r-^^-\ (10) 

where the corresponding excitation is given by (1;0), and consequently by Pb,x = A/2. 
Due to the selection rule (^) for fermions, the corresponding excitation of the correlation 

(0U^)0f(O))a oc cos(A;^r) r"^^-^ (11) 

is given as (1, |) and thus Pf,x = (A + l/A)/2. 

Taking the Fourier transform of the long-distance expressions (|10D and (pH]) , we get 
the momentum distribution function n{k). The long-distance behavior now translates into 
the behavior at long-wavelength or small wave numbers, and thus for the Bose momentum 
distribution near the origin, 

nsik) oc = \kf^''-\ (12) 

For the Fermi case, we have correspondingly for the momentum distribution near kp, 

npik) (x\k- A;i.|"^-^sgn(fc - kp) = \k - kpf^'^'^sgn^k - kp) (13) 

In Fig. |1] we show the behavior of the critical exponents as functions of A. We emphasize 
that the exponents do not depend on the particle density d, but on A only. 

Our above derivation of the critical exponents follows largely the article of Kawakami 



and Yang [jl4[. Other derivations have been given by Mironov and Zabrodin, who obtained 
the correlation exponents for a wider class of ID quantum models with pairwise interaction 
||18|| , and by Kolomeisky |jl5[, who applied methods of two-dimensional elasticity theory to 



arrive at corresponding results. 



III. CORRELATIONS IN THE SUTHERLAND MODEL 

The ground state wave function for the periodic g/r'^ model with M particles, either 
bosons or fermions, interacting on a one-dimensional ring of circumference of length L, is of 



the Jastrow product form and given by p2,23 



= Cx^M n V'A [27r(xfe - Xj)/L] . 

l<j<k<M 

The factor 5* is equal to 1 for bosons, while for fermions S = (—1)^ is the parity of the 
particle ordering permutation P, xpi < ■ ■ ■ < xpM- The periodicity in L is assured if M 
is odd. ipxiO) = |2sin(^/2)|'^ for bosons and ip\{0) = |2sin(^/2)|''^sgn(6') for fermions, the 
normalization constant is given as 
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r2 _ r(i + A)^ 



As has been noted in |2^, the square of ^\ for the "special" values g = — 1,0,4 — or 
A = |, 1, 2 — may be recognized as being identical to the joint probability density function for 
the eigenvalues of matrices from Dyson's circular ensemble [Q, A = |, 1 and 2 corresponding 
to orthogonal, unitary and symplectic ensembles, respectively. 

Results from the theory of random matrices now enable the calculation of various cor- 
relation functions [0. We will use the following result: For u{9) = u{—9) = u{2ti + 6), 
arbitrary otherwise, and the special values of A, 

M 

(n«(^.))A = det[F,,], (14) 



where for A = 1, 



Fpi = ^ r ^H^) cos(p - q)9, 



p,q = |,|,...,M- |; 



for A = 2, 



1 

Fpg = / d9u{9) [pcos{p9) cos{q9) + qsm{p9) sin(g^^)] 



2ttp J-tt 



' ' M-h 



finally for A = |, 
- / 



Fpq = -T- j d9u{9) I d(f)u{(j))u{(f))sgn{9 — 0) [cos{p(f)) sm{q9) — cos{p9) sin(g0)] 



i^, y 2;2'''''2 2 

and we assume M to be even for all three cases. Thus, Fpg is an M x M matrix for A = 1, 2 
and an Y ^ Y i^^-trix for A = |. 

The one-particle density matrix is defined as 

M-l L 




M-l L 

p{x-x') = MYI / dxj-^x{xi,...,XM-i,x)-^x{xu...,XM-i,x'), 



which with the help of (|1^) may be written as 



(J2 M-l 

p{x - x') = M-^{ n ^a(% - 2nx/L)i,x{9, - 2t,x' / L))x,m-i 

^A,Af-l j=l 

_ M r(i + A)r(i + A(M-i)) 

" L r(i + AM) '^''^^^^^J- ^^^^ 
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In the last formula, we use the appropriate matrix Fpq with M replaced by M — 1 . By the 
periodicity condition on the wave function, we restrict M to be odd. Thus, M — 1 is even 
and our above formulae are consistent. The normalization of p has been chosen such that 
p(0) = M/L = d. The momentum distribution n{k) is defined as the Fourier transform of 
p, i.e. 



n{k) = / drp(r)e ''''' 

J —oo 

and the normalization is such that 

/oo 
dkn{k). 
-oo 

The A = 1 case corresponds to free fermions or hard core bosons. Free fermions are well 
understood, and we find 



^(,) = M-dr)^ ^^^^ ^ \ ^, \k\ < nd 
nr 



27r' 

0, \k\ > ird 



Thus, ap,! = in agreement with (13). Hard core bosons have been treated extensively 
in the literature |jTl],^,|T6| . There is no condensation into a single momentum state, and 
it is proven that n{k) diverges as \k\~2 at the origin. Again, this is in agreement with the 
corresponding result of section ||. 

For A = 2, the matrix Fpq has a particularly simple Toeplitz form, which allows the 



explicit calculation of p p2 



p(r) = "P""->. Si(r)^ rd.?'°W 



27rr Jo X 



n(k) = 4. ' f - (16) 

Thus, the conformal result aB,2 = agrees again; although, the logarithmic behavior is not 
universal and so does not follow from the finite-size analysis. 

This leaves the remaining three cases — bosons with A = |, and fermions with A = | 
and A = 2 — to be investigated numerically. 

Let us now study the restriction of the above family of wave functions to a lattice. The 
size of the systems L is now an integer iV, and the coordinates of the particles are restricted 
to integers Xj = 1,2, . . . , N. The wave functions are nearly as before; the only substantial 
change is the normalization constants Ca,m- For the special values of A, the new C^m 



are known |2^. The one-particle density matrix p(r) and the corresponding momentum 
distribution function n{k) for the lattice case may still be calculated as before, if we replace 
the integrals by appropriate sums over lattice sites. 

Note, however, that whereas the density d = M/ L in the continuum case entered trivially 
into all expressions as a length scale, the new density d' = M/N in the discrete case is 
dimensionless and will thus enter all expressions in an essential way. Indeed, the previous 
continuum results are the zero-density limit of the discrete case. We showed in section |I| that 
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the critical exponents of the momentum distribution, as obtained from conformal analysis, 
should be independent of density. This prediction may therefore be verified by studying 
both continuum and lattice cases. 

However, we want to emphasize that whereas the wave functions for the continuous case 
describe the ground state of the periodic g/r"^ model, there need not be any corresponding 
Hamiltonian for the wave functions of the lattice gas. For A = 1, we of course have either 
hard core bosons of free fermions; while for A = 2 bosons, we have the Haldane-Shastry 
model of a Heisenberg magnet with exchange. However, for the remaining 



three cases — precisely the cases we are to evaluate numerically — it is not clear to what 
specific model our results belong to in the discrete case. From the conformal viewpoint, of 
course, the results represent a universality class of models. 

Finally, the lattice gas can be shown to exhibit a particle-hole symmetry, such that we 
only need to consider d <^. In addition, at half- filling, Fermi statistics ensure p(r) = for 
r 7^ or even and therefore for the momentum distribution, the symmetry n{k)+n{k+Ti) = 1 
holds. 



IV. RESULTS 

Using the results obtained in section |T|, we calculate the one-particle density matrix 
p{x — x'), or by translation invariance p(r), with the help of (^) in terms of the appropriate 
matrix Fpg. Taking the Fourier transform, we compute the momentum distribution n{k), 
which, for fixed L and N, turns out to have a slightly different quantitative behavior for the 
continuum and discrete cases. In Fig. Q and |^, we show n{k) corresponding to the discrete 
case at the special values A = |, 1, 2 and half-filhng, for bosons and fermions, respectively. 

The extraction of the critical exponents a^.A for Bose and ap^x for Fermi statistics is 
complicated by finite-size effects. The interchange of the L — > cx) and A; — > (/c — > kp) 
limits is non-trivial, and the correct order — L — oo limit first — is hard to achieve. In 
addition, we know from ([T6|) that the momentum distribution for the bosonic A = 2 case has 
a logarithmic behavior, which is not predicted by the conformal analysis, and the possible 
appearance of this or some other finite-size effect for the remaining cases may further obscure 
the extraction of the correlation exponents. 

However, a careful analysis of our data reveals that the one-particle density matrix p 
obeys an unusual scaling law. For different lattice sizes L, a limiting curve poo may be 
constructed as 

Poo(r) = ^PLir/L), (17) 

both for Bose and Fermi statistics. Therefore, the limit L ^ oo of pL(r) may be taken at 
any point r, and we can directly extract the asymptotic value of the critical exponents of 
the momentum distribution by studying the long-distance behavior, i.e. p(L/2). In Fig. ^ 
we show the remarkably good convergence of our estimates for the critical exponents as 
function of the inverse length 1/L of the ring for the two A = | cases and the fermionic 
A = 2 case. 

The scaled one-particle density function exhibits another quite remarkable property. In 
Fig.'s ^ and ^, we have plotted poo for bosons and fermions and normalized it so that its 
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value at the end points equals 1. Although, as we remarked above, n{k) is different for the 
discrete and the continuous case, poo is identical for both cases up to the accuracy of our 
calculation and thus independent of density. In addition, the curves corresponding to A = ^ 
and A = 2 in Fig. ^ lie on top of each other. This seems to indicate that poo is the same 
for identical values of the critical exponent a. This is further supported by the fact that 
the bosonic A = 2 curve, corresponding to a = 0, is the same as the fermionic A = 1 case, 
corresponding to a = 0, too. We will discuss the implications of these findings in more 
detail in the next section. 

We have calculated the critical exponents, both for the bosonic and fermionic cases and 
A = 1, 2, for up to M = 641 particles on L = 1282 sites. Since the scaling curves are identical 
for discrete and continuum cases, so are the exponents. An independent estimate from n{k) 
directly confirms this. For A = |, the continuum formulae require the computation of double 
integrals of transcendental functions rather than orthogonal functions, and so do not reduce 
to the calculation of simple Kronecker 5's. Therefore, calculations could be done only up 
to M = 41 particles at half-filling. The discrete expressions require less computational 
effort and data for up to M = 161 for the bosonic case, and M = 321 for the fermion case 
have been calculated. Again, the results indicate that the critical exponents are the same, 
independent of whether the continuum or the discrete formulae are used. Since the discrete 
case is the low-density limit of the continuous case, we have therefore established that up to 
our numerical accuracy, the critical exponents aF,x and aB,x are independent of the density 
d as predicted by the conformal arguments. 

Let us now compare the values of our numerical estimates with the predictions of section 
^ for the critical exponents. Preliminary estimates for small particle numbers indicated 
a disagreement with the predictions of the conformal analysis for the bosonic A = | case 
2B| . In table |, we show that our present study shows excellent agreement of conformal and 



numerical results; the estimates for the exactly known results are included as a check on the 
numerical analysis. 

Finally, we note that we could not find any logarithmic behavior except — as predicted 
— for the bosonic A = 2 case. However, there do seem to be substantial finite-size effects. 



V. DISCUSSION 

In the preceding sections, we have reviewed arguments of conformal finite-size scaling 
analysis leading to predictions for the critical exponents of the momentum distribution in 
the ID quantum liquid with long-range interaction g/r"^. Both for bosonic and fermionic 
systems and at the special values g = —|, 0, 4 for the interaction strength, we then proceeded 
to calculate the same exponents by use of explicit formulae as given by random matrix 
theory. For non-trivial cases, no analytic solutions exist and we calculated the exponents 
numerically. This is most conveniently done by discretizing the continuous problem to a 
lattice problem. The continuum results are then recovered by taking the low-density limit 
of the lattice system. Our results show that the values of the correlation exponents are 
the same for lattice and continuum case, thus independent of density. This verifies the 
corresponding prediction of conformal field theory. In addition, the numerical estimates of 
aB,<yF are in excellent agreement with the conformal predictions. 
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The numerical analysis was helped by the appearance of the unusual scaling law (p!7D. 
As it turns out, the scaled and properly normalized one-particle density function may 
be parametrized by the critical exponent a only, independently of statistics and interaction 
strength A. Furthermore, the interacting Bose gas at A = 2 has the same scaled poo as 
the non-interacting free Fermi gas. This seems to suggest that not only the strength of 
the interaction but also the specific type of the interaction is irrelevant. We do not know 
whether this behavior may be explained in terms of the conformal symmetry of the system, 
but hope to come back to this point in a future publication. 

Trying to classify the g/r'^ model by its critical anomaly, we calculated the finite-size 
corrections to the ground state energy as in (|1|). We found that in the cylindrical geometry, 
c depends continuously on the value of the coupling constant A. On the other hand, use 
of the low-temperature expansion of the free energy showed that c = 1. Furthermore, the 
exponents a^, ap vary continuously as functions of A. In [H this has been shown to be a 
typical phenomena of conformal field theories with anomaly c = 1. In addition, Kawakami 
and Yang [|1^] demonstrated that [3b, Pf in our model obey the scaling relations of the c = 1 



Luttinger liquid. Therefore, we believe that c = 1 is the correct value of the critical anomaly 
of the g/r'^ model. We emphasize that for models with c > 1 there is no Kac- formula, 
relating the critical anomaly and the scaling dimensions, so we proceeded to calculate the 
correlation exponents independently of c. 

We conclude this paper by noting that at present, we do not know a satisfactory rigorous 
explanation for the observed discrepancy in the estimate of c in terms of conformal field 
theory. However, we believe this to be a remainder of the long-range character of the g/r"^ 
interaction, and in the light of our introductory remarks of section |, it should not be too 
surprising if time correlations differ from space correlations in long-ranged models. 
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FIGURES 



FIG. 1. Correlation exponents for bosonic and fermionic systems as functions of the interaction 
strength g. The values corresponding to the special points g = —^,0,4 are indicated. 

FIG. 2. Momentum distributions for the bosonic system on a lattice at half-filling, for the 
special values of A. 

FIG. 3. Momentum distributions for the fermionic system on a lattice at half-filling, for the 
special values of A. 

FIG. 4. Estimates of the critical exponents as functions of the inverse length of the lattice. (□) 
corresponds to /3d i, (+) to 0:^2 and (O) to a^i. 

'2 ' '2 

FIG. 5. Scaled one-particle density poo for the bosonic system at half-filling, for the special 
values of A. The curves for the discrete and the continuous case are identical up to the accuracy 
of the plot. The end points are normalized to 1. 

FIG. 6. Scaled one-particle density poo for the fermionic system at half- filling, for the special 
values of A. The curves for the discrete and the continuous case are identical up to the accuracy 
of the plot. The end points are normalized to 1. Note that the curves for A = ^ and A = 2 are 
identical, too. 



12 



TABLES 



TABLE I. Results for the critical exponents a from conformal finite-size scaling analysis and 
numerical calculations for the special values of A. The estimated error for the last digit is given in 
parenthesis. The * indicates that the value of the corresponding exponent is known exactly. 
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